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[57] ABSTRACT 

The invention is a ride quality meter that automatically 
transforms vibration and noise measurements into a 
single number index of passenger discomfort. The noise 
measurements are converted into a noise discomfort 
value. The vibrations are converted into single axis 
discomfort values which are then converted into a com- 
bined axis discomfort value. The combined axis discom- 
fort value is corrected for time duration and then 
summed with the noise discomfort value to obtain a 
total discomfort value. 

8 Claims, 7 Drawing Figures 
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RIDE QUALITY METER 

ORIGIN OF THE INVENTION 

The invention described herein was made by employ- ^ 
ees of the U.S. Government and may be manufactured 
and used by or for the Government for governmental 
purposes without the payment of any royalties thereon 
or therefor. 
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BACKGROUND OF THE INVENTION 

The invention relates generally to passenger ride 
discomfort and more specifically concerns a tool for 
estimating passenger ride discomfort within complex 
ride environments. 15 

The prior art methods and devices for measuring 
passenger ride discomfort generally measure and dis- 
play individual vertical or lateral accelerations that are 
frequency weighted according to ride comfort curves 
derived from motion sickness data and the International 20 
Standards Organization recommendations. These de- 
vices and methods output a weighted value of accelera- 
tion for either of the two axes of motion. The resultant 
weighted acceleration level is usually interpreted in a 
dichotomous manner, i.e., representing a ride that is 25 
comfortable or uncomfortable. 

The disadvantages of the prior art include the follow- 
ing. The output of existing devices and methods is a 
frequency weighted value of physical acceleration 
which does not necessarily relate on a one-to-one basis 30 
with the associated subjective experience of the mea- 
sured environment. For example, many different vibra- 
tions (which produce large variations in subjective dis- 
comfort) can have identical weighted acceleration lev- 
els. Thus, a weighted output of these devices and meth- 35 
ods does not generally represent a single unique level of 
subjective discomfort/acceptance. Existing devices and 
methods are limited to the production of comfort indi- 
ces for single axes of vibration. They cannot provide 
reliable estimates of discomfort due to simultaneous 40 
vibrations in more than one axis of vibration. They 
canot adequately assess the effect upon subjective dis- 
comfort of multiple frequencies of vibration within 
individual axes; they do not account for the interaction 
and/or summation of the effects of combined noise and 45 
fibration upon human discomfort/acceptance; they do 
not incorporate the effects of angular vibrations, either 
singly or in combination with linear vibrations; and they 
do not account for adaptation of passengers to the vi- 
bration environment. 50 

It is an object of this invention to provide a ride qual- 
ity meter in which simultaneous measurements of inte- 
rior noise and vibrations in the transportation system 
are used. 

Another object of this invention is to utilize the com- 55 
bined noise and vibration measurements, together with 
empirically derived psychophysical laws governing 
human discomfort response to combined noise and vi- 
bration to generate and display, in real time, an index of 
passenger discomfort that directly relates to passenger 60 
subjective acceptance of the measured ride environ- 
ment. 

A further object of this invention is to display at the 
option of the user, the contribution of the various com- 
ponents of the ride environment to the total discomfort 65 
experienced by passengers. 

Still another object of this invention is to provide a 
ride quality meter in which the effects of multiple fre- 
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quencies and multiple axes of vibration are automati- 
cally accounted for within the meter. 

Other objects of this invention will become apparent 
hereinafter in the specification and drawings. 

SUMMARY OF THE INVENTION 

the invention is a ride quality meter that provides a 
value which is indicative of the discomfort that a pas- 
senger will experience during a ride on a passenger 
vehicle. 

A transducer box containing accelerometers is placed 
on the vehicle such that electrical signals are produced 
that are proprotional to the vehicle vertical, lateral, 
longitudinal, roll and pitch accelerations. A micro- 
phone on the vehicle produces an electrical signal pro- 
portional to noise. The electrical signal proportional to 
noise is divided into several octave bands and then 
applied to a noise discomfort computer which computes 
a noise discomfort value. 

The electrical signals from the accelerometers are 
applied through a fast Fourier transform signal analyzer 
and computers to obtain single axis discomfort values. 
These values are combined by computer means to ob- 
tain a combined axis discomfort value. The combined 
axis discomfort value is corrected for time duration and 
then summed with the noise discomfort value to obtain 
a total discomfort value which is displayed on a display 
means, 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a graph of discomfort as a function of per- 
cent of passengers uncomfortable; 

FIG. 2 is a block diagram of the concept of the inven- 
tion; 

FIG. 3 is a graph showing the time duration correc- 
tion needed to be made to the combined axis discomfort 
value; 

FIG. 4 is a graph of noise corrections as a function of 
the discomfort level due to the fibration component of 
the ride environment; 

FIG. 5 is a schematic drawing of the use of the in ven- 
tion on a passenger vehicle; and 

FIG. 6 (a) and (b) is a block diagram of the invention. 

DETAILED DESCRIPTION OF THE 
INVENTION 

This invention is the result of a program utilizing 
approximately 2200 test subjects in the development of 
a generalized ride comfort model for engineering de- 
sign/analysis applications. The model provides a single 
numerical descriptor of passenger discomfort which is 
measured along a ratio scale of discomfort that is an- 
chored at discomfort threshold. This discomfort scale is 
illustrated in FIG. 1, which shows the relationship be- 
tween the discomfort scale (ordinate) and the core- 
sponding percentage (abscissa) of passengers who 
would rate that discomfort level as being uncomfort- 
able. A value of unity along the discomfort scale corre- 
sponds to discomfort threshold, i.e., 50% of the passen- 
gers would be uncomfortable. 

The ride comfort model concept of this invention is 
illustrated in FIG. 2. The model 7 shown enclosed by 
the dashed lines consists of three basic elements; (1) 
empirical estimation of discomfort due to sinusoidal 
and/or random vibrations within single axes 8; (2) em- 
pirical estimation of the discomfort due to vibration in 
combined axes 9; and (3) application of empirically 
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determined corrections for the effects of interior noise 
and duration of vibration 10. Input to the model is the 
vehicle ride environment, and output of the model is the 
total discomfort measured along the discomfort scale of 
FIG. 1, Each of the model elements is discussed in detail 
below. 

Single axis discomfort to sinusoidal vibration: Dis- 
comfort responses of passenger subjects to vibrations 
applied in each of five axes of motion (vertical, lateral, 
longitudinal, roll, and pitch) were obtained and mod- 
eled in a recent study. Responses to sinusoidal vibration 
were obtained only for the vertical, lateral, and roll axes 
of motion, since sinusoidal pitch and longitudinal vibra- 
tions were not considered typical of actual vehicle ride 
environments. The equations developed for describing 
discomfort response to sinusoidal motions are given by: 

Ds=-a + bXpJovXp^X (la) 


4 

TABLE 1-continued 


Frequency, Hz 

a 

b 

c 

8 

-0.4184 

14.8952 

7.92 

9 

-0.0636 

11.6969 

10.64 

10 

0.3307 

8.9291 

14.44 


(c) Roll Sinusoidal 


1 

-2.31 

5.85 

1.239 

2 

-0.18 

4.70 

N/A 

3 

0.28 

2.50 

N/A 

4 

0.35 

2.35 

N/A 


Single axis discomfort to random vibration: It is conve- 
nient to express the discomfort response to random 
vibration in the vertical and lateral axes in the form of 
15 single equations obtained from multiple regression anal- 
yses. This provides a means for estimating discomfort 
response using bandwidths and center frequencies. The 
single relationships for random vertical vibration are 


Ds=cXp, Xp<X 


(lb) 20 D>r/?=-1.75+O.857(Ci^-O.lO2(C/02+O.OO346(C- 

F)^ + ^^A%rms. m 


where 

D^— discomfort due to a sinusoidally applied vibra- 
tion; 

X^=peak linear (in g units) or angular (rad/s^) for 25 
roll acceleration level; 

X=0.06 g for linear (vertical and lateral) sinusoidal 
acceleration, 0,50 rad/s^ for roll acceleration; and 

a,b,c,= empirical constants which depend upon the 
axis and frequency of the sinusoidal vibration. 30 
Values of a, b, and c, the coefficients of equations (la) 
and (lb), are given in Table 1 below for sinusoidal verti- 
cal, lateral, and roll vibration, respectively. 


TABLE 1 


Frequency, Hz 

a 

b 

c 

1 

(a) Vertical Sinusoidal 

.3946 8.8296 

15.41 

2 

-.3713 

15.2731 

9.08 

3 

-.7685 

21.4441 

8.64 

4 

-1.0028 

27.1273 

10.41 

5 

-1.2352 

32.2146 

11.63 

6 

-.7592 

28.8279 

16.17 

7 

-.7188 

27.4856 

15.51 

8 

-.0576 

19.8988 

18.94 

9 

-.8919 

21.9987 

7.31 

10 

-1.2718 

22.9530 

1.76 

11 

-.6912 

16.9931 

5.47 

12 

-.4937 

14.0437 

5.82 

13 

-.3695 

12.0297 

5.87 

14 

. -.3470 

10.7501 

4.97 

15 

-.5220 

10.4234 

1.72 

16 ‘ 

-.1406 

8.3656 

6.02 

17 ' 

.1650 

6.899.7 

9.65 

18 

-.2190 

7.5948 

3.94 

19 

-.3326 

7.5326 

1.99 

20 

.0986 

6.1421 

7.79 

21 

-.1989 

6.7045 

3.39 

22 

-.1769 

6.5021 

3.55 

23 

.0345 

5.9102 

6.49 

24 

-.0465 

6.0773 

5.30 

25 

.0494 

5.8456 

6.67 

26 

.0010 

6.0208 

6.04 

27 

-.0684 

6.2664 

5.13 

28 

-.1695 

6.6472 

3.82 

29 

-.0324 

6.4483 

5.91 

30 

-.0766 

6.7358 

5.46 

1 

(b) Lateral Sinusoidal 
-0.8322 26.7849 

12.91 

2 

-1.1106 

52.2679 

33.76 

3 

-0.3586 

32.1940 

26.22 

4 , 

0.0217 

19.9130 

20.27 

5 

-0.3163 

19.0267 

13.76 

6 

-0.7048 

19.8629 

8.12 

7 

-0.7024 

16.3704 

4.66 


35 


40 


45 


50 


55 


60 


65 


and for random lateral vibration 

D^'/?=0.89-0.1S7(Cf)+0.016(^ffO + 29.15 grms (3) 

where 

DffT?, discomfort response to fandom vertical, 

lateral vibration; 

CF= center frequency in Hz of applied random vi- 
bration; 

^W=bandwidth in Hz of applied random vibration 
(defined by 10-dB downpoints); and 

gr/ 725 = root-mean-square level of applied vertical or 
lateral random vibration with each defined band- 
width 

Equation (2) is valid for center frequencies in the range 
of 2 to 13 Hz and bandwidths over the range of 2 to 10 
Hz whereas equation (3) applies to center frequencies in 
the range of 2 to 9 Hz and bandwidths of 2 to 10 Hz. 

Discomfort functions for random roll, pitch, and 
longitudinal vibrations are given by the following: 

DBR=^0M~^4M &rms\ for &V/«5=0-141 rad/s^ (4) 

£>0/?=7.O4 Orms'.foT 0.141 rad/s^ 

/?4>/?=0.4I+5.07 (frms; for $rms^0.U6 rad/s^ (5) 

Dff,R^8.62 $rmsi for $rm5<0.tl6 rad/s2 

DuR^ -0.02+42.24 (g,,„,)u (6) 


where 

Da/?, Dtf)/?, Du/? = the discomfort due to random vi- 
bration in the roll, pitch, and longitudinal axes, 
respectively; 

root-mean-square acceleration level, in rad/s^, 
of a random roll vibration having a bandwidth of 5 
Hz and centered at 3 Hz; 

root-mean-square acceleration level, in rad/s^ 
of a random pitch vibration having a bandwidth of 
5 Hz and centered at 3 Hz; and 
(gr/w 5 )w= root-mean-square acceleration level, in g, of 
a random longitudinal (±g) vibration having a 
bandwidth of 5 (or 10) Hz and centered at 5 Hz. 
Combined-axes discomfort: Since many transporta- 
tion vehicles contain vibrations in more than one axis at 
a time, a series of experiments was conducted to de- 
velop a reasonable procedure for estimating discomfort 
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response to the combined-axes situation. Subjective 
reactions to combined vertical, lateral, roll vibrations 
were obtained from 126 subjects, and 54 subjects were 
used to obtain subjective reaction to combined vertical, 
longitudinal, pitch vibrations. The reason for the em- 5 
phasis (in terms of subjects) oh the vertical, lateral, Toll 
combination arose from the fact that it is the most im- 
portant in terms of vehicle ride quality. 

The characteristics of the physical stimuli used in the 
combined axes experiments are given in Table 2 below lo 
for both axis combinations. 


example, a model based upon physical descriptors such 
as vertical, lateral, and roll acceleration level would 
permit adequate estimation of discomfort response only 
for the particular vibration frequency spectra that were 
used to derive the model. Application of such a model 
to vibratiohT haying different frequency characteristics 
could produce large errors iri estimated discomfort. 

Duration of vibration: An extensive study using 210 
passenger subjects determined correction factors that 
were incorporated into the ride comfort modeling pro- 
cess to account for the effects of vibration duration. 


TABLE 2 


Axis 

Com- 

bina- 

tion 

Range of Physical Stimulus Values 
Used in Combined Axis Studies 

Bahdwidth, 

Hz 

Axes 

Included 

Range of 

Acceleration 

Level 

Center 

Frequency, 

Hz 

1 

Vertical 

0.025-0.075 g 

3,5,7 

2,5,10 


Lateral 

0.025-0.075 g 

3,5,7 

2,5,10 


Roll 

0.02-0.50 rad/s2 

3 

5 

2 

Vertical 

0.025-0.075 g 

5 

10 


Longitudinal 

0.025-0.075 g 

5 

10 


Pitch 

0.02-0.50 rad/s2 

3 

5 


Their results indicated that subjects adapted to vibra^ 
tion ride environments haying durations up to one hour. 
This adaptation process was independent of acceleratin 
level and varied linearly with vibration duration, i.e., 
perceived discomfort decreased linearly with increasing 
yibratiop exposure ^ time. . The i:elaiionships describing 
the decrease in discomfort as a function of time for 
durations of up to 120 minutes are given by; 

/ ' A£)£,= d 003 1 - 0.0 1 2/; for 1 60 minutes (i 3) 

and " r ' , / ' ^ / r . 


The subjective responses, together with the measured 25 
values of the physical stimulus factors, were used as 
input to both linear and polynomial multiple regression 
routines that computed various least-squares models to 
fit the empirical data. The model selected as best for 
general use in estimating discomfort response to various 30 
combinations of vertical, lateral, and roll vibrations is 
given by; 

-0.44+ 1.65 />c/;forT?c/^0.88/? (7) 

35 

for Dc/<0.88 D (8) 

where 

Da=y/b^w+D^v+D^e ( 9)40 

Where Dff', Y>Vy D^ are the discomfort levels due to 
multiple frequencies (computed by equation (19) below) 
within the vertical, lateral, and roll axes, respectively. 
Discomfort due to combined vertical, longitudinal, and 
pitch axes of motion can be estimated by the following 
equations; 

£>W^u({,= -1.07+1.77 Dc2;forDc2^1.0i7 (10) 

£>H'u*= 0.70 Dcj: for Dc 2 < 1.0 0 (11)50 

where ; 

Dci = ^Dn^ + D^„R + b^i,R (12) 

55 

Equations (7) and (10) are based upon experimental data 
which produced values of Dc/ in the range of 0.88 to 
5.00 D and values of Dc2 over the range 1.00 to 5.00. 
Equations (8) and (1 1) represent reasonable approxima- 
tions of discomfort response below discomfort thresh- 60 
old. The advantage of modeling combined axis discom- 
fort response in terms of discomfort due to each individ- 
ual axis of vibration is that computation of individual 
axis discomfort inherently accounts for the effects of 
vibration frequency. Fn addition, it increases the genef- 65 
ality of the model, since values of D>f 7 ?, D^t?, etc., 
represent subjective units of discomfort which can re- 
sult from any number of different ride spectra. For 


AZ>xj=— 0:72; for 60 minutes 120 minutes (14) 

'Where- -'"'" ^ ■ ■ ' 

ADi?=the duration correction at time t, and 

t=the duration in minutes of vibration. 

The vibration duration correction corresponding to 
equations (13) and (14) is shown in FIG. 3. Note that at 
t= 60. minutes, the vibration duration correction is ap- 
proximately —0.72 D. The corrections given by equa- 
tions (13) and (14) provide reasonable estimates of the 
duration corrections for trip times up to 120 minutes. 
The assumption of a constant duration effect over the 
range of 60 to 120 minutes was believed to be justified 
on the basis of the fact that no significant effect of fa- 
tigue or further adaptation has been observed* in any 
ride quality experiment that utilized subjects for test 
periods of two hours or more. Using the above correc- 
tion, the total vibration discomfort, corrected for the 
effectofduration,.isgivenby; 

Dvib-Dc+t^Do (15) 

where 

D^ 7 (&= total vibration discomfort in D corrected for 
the effect of duration; and 

Dc^ vibration discomfort in D uncorrected for the 
effect of duration. ^ , 

Combined noise/vibration: The approach used to 
define discomfort due to a combined noise and vibration 
environment was to expose passenger-subjects to vari- 
ous parametric combinations of noise and vibration and 
to obtain their subjective reactions to the total environ- 
ment. These studies indicated that the discomfort due to 
the combined* noise and vibration environment could be 
represented as a summation of discomfort components 
attributable to each of the physical components. Total 
discomfort in the combined environment was given by 
the following: 

DTotal=^Vib+^N (16) 

where 

D 7 b/o/=the discomfort, in D, due to the combined 
noise and vibration environment; and 
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DAT=the contribution to total discomfort attributable 
to the noise, in D, present within a vibration envi- 
ronment. 

The following paragraphs discuss the equations and 
procedures developed to estimate the noise discomfort 5 
contribution when the noise spectrum contains either a 
single octave band or contiguous octave bands. 

Noise discomfort contribution — single octave bands: 

A recent study derived a model of subjective discom- 
fort response to combined noise and vibration that en- 10 
compassed six octave bands of noise (63 to 2000 Hz) and 
a wide range of sinusoidal vibration. The metric used in 
that investigation was the A-weighted sound pressure 
level. It was found that noise discomfort varied with 
noise level, noise octave frequency, and the level of 15 
vibration present in the environment. Thus, it was nec- 
essary to model the interactive effects of the two physi- 
cal kimuli in order to provide accurate estimation of 
total discomfort in the combined environment. The 
interactive effects of noise and vibration are illustrated 20 
in FIG. 4 which shows the noise discomfort correction 
as a function of vibration discomfort level for several 
noise levels. The noise discomfort correction varies in 
an approximately linear fashion with the level of vibra- 
tion discomfort present in the ride environment. Conse- 25 
quently, the model selected to represent this effect and 
to estimate noise corrections due to a single octave band 
of noises is given by: 

DN{jj)=ai-\-biDyibWFi ( 17 ) 39 

where 

DAr(o)= noise discomfort due to the ith octave band 
having an A-weighted noise level of j dB in the 
presence of a vibration environment that produces 
a vibration discomfort level of Dvw 
a/,b/,= empirically determined coefficients; and 
WF/=a weighting factor that corrects for the effect 
of the ith noise octave band. The weighing factors 
are given in Table 3. 

TABLE 3 

Octave Band Weighting Factors 


Octave Center Frequency, Hz Weighting Factor 


63 

1.470 

125 

0.963 

250 

0.786 

500 

0.646 

1000 

0.688 

2000 

1.448 


Values of a/ and b/are given in Table 4 for A-weighted 
noise levels ranging from 65 to 100 dB(A) and vibration 
discomfort levels from 0.5 to 4.0 D. 

TABLE 4 


Values of Slope and Intercept for Equation 17 


Noise 

Level 

dB(A) 

Intercept 

a 

Slope 

b 

Noise 

Level 

dB(A) 

Intercept 

a 

Slope 

b 

65 

0.3447 

-0.1219 

83 

2.2294 

-0.5118 

66 

0.4172 

-0.1445 

84 

2.3718 

-0.5329 

67 

0.4935 

-0.1669 

85 

2.5164 

-0.5533 

68 

0.5736 

-0.1893 

86 

2.6649 

-0.5738 

69 

0.6575 

-0.2116 

87 

2.8172 

-0.6145 

70 

0.7452 

-0.2337 

88 

2.9732 

-0.5942 

71 

0.8368 

-0.2558 

89 

3.1330 

-0.6346 

72 

0.9320 

-0.2777 

90 

3.2968 

-0.6547 

73 

1.0312 

-0.2995 

91 

3.4642 

-0.6746 

74 

1.1340 

-0.3212 

92 

3.6354 

-0.6944 

75 

1.2408 

-0.3429 

93 

3.8104 

-0.7142 

76 

1.3512 

-0.3644 

94 

3.9893 

-0.7338 
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TABLE 4-continued 


Noise 

Level 

dB(A) 

Values of Slope and Intercept for Equation 17 

Intercept 

a 

Slope 

b 

Noise 

Level 

dB(A) 

Intercept 

a 

Slope 

b 

77 

1.4654 

-0.3858 

95 

4.1720 

-0.7533 

78 

1.5835 

-0.4071 

96 

4.3574 

-0.7724 

79 

1.7055 

-0.4284 

97 

4.5486 

-0.7921 

80 

1.8311 

-0.4494 

98 

4.7426 

-0.8113 

81 

1.9605 

-0.4704 

99 

4.9404 

-0,8304 

82 

2.0938 

-0.4913 

100 

5.1421 

-0,8494 


Note that the noise discomfort contribution to total 
discomfort response requires that the discomfort due to 
the vibration components of the ride environment be 
computed first. This is done with the use of the equa- 
tions and procedures discussed earlier. 

Noise discomfort contribution — contiguous octave 
bands. The relationship used to compute the noise dis- 
comfort contribution to the total discomfort response 
when noise is present is more than one octave is given 
below: 

Z)yV “ +0.3 (18) 

where 

Dyv=the noise discomfort due to a continuous noise 
spectrum in the presence of a vibration environ- 
ment that produces a discomfort level of Y>vib- 
Note that values are computed from equation 

(17) and Div(y-)^a;c represents the noise octave band that 
provides the greatest discomfort. 

Equation (17) or (18) is sufficient to compute an esti- 
mate of the noise discomfort contribution to the total 
discomfort response. Application of equation (16) is 
then sufficient to estimate total discomfort response, 
and FIG. 1 can be used to determine the percentage of 
passengers finding that particular ride environment 

, uncomfortable. 

40 

Multiple frequency vibrations: An accepted method 
for handling a vibration spectrum containing multiple 
frequency components is not presently available. The 
International Standards Organization (ISO, 1972) rec- 
ommends that, for both discrete and narrowband vibra- 
tion, the rms acceleration within each band (or fre- 
quency) be evaluated with respect to the appropriate 
limit of that band (or frequency). The implicit assump- 
tion in the ISO approach is that acceptability of a given 
jQ ride environment is determined by the dominant com- 
ponent of the vibration spectrum. This assumes that no 
interactions occur between the discomfort produced by 
different frequencies. A recent study directly addressed 
the effects of multiple frequency vertical (±gz) vibra- 
55 tion upon subjective discomfort and examined the ap- 
propriateness of several candidate models for depicting 
such effects. Using a discomfort matching procedure, it 
was determined that the discomfort of many multiple 
frequency motions may be meaningfully expressed in 
60 terms of the equivalent level of a single frequency vibra- 
tion. Specifically, the levels of a 10 Hz vibration that 
were equivalent to complex vibrations were found to be 
well predicted by the rms levels of the 10 Hz equivalent 
to the individual sinusoids present in the complex vibra- 
65 tion. This approach is directly analogous to a vector 
summation of discomfort units due to individual fre- 
quency components within a single axis ride spectrum. 
It is, therefore, assumed that adequate prediction of 
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discomfort due to multiple frequency components can 
be obtained from the following equation. 



where 

D»^^=the within axis discomfort due to one or more 
narrowband random and discrete frequency com- 
ponents; and 

D/=the discomfort due to the ith narrowband or 


10 

square level {gwRy gVRr gBRyZuRy gi^R) of each compo- 
nent. The spectrum identification module is a software 
element that scans the five spectral density digital out- 
puts and performs the following functions: (1) identify 
the single maximum (peak) level within each spectrum; 
(2) estimate the bandwidth of the largest spectrum peak 
by defining it as the frequency range contained within 
the 10 dB downpoints of the peak; (3) identfy the center 
frequency as the midpoint frequency within the defined 
bandwidth; and (4) compute the root-mean-square ac- 
celeration level contained within the defined bandwidth 


discrete frequency component within a given axis 

(computed from equations (1) to (6), where appro- 

priate). . . is 

Turning now to the embodiment of the invention 
selected for illustration the number 11 in FIG. 5 desig- 
nates a vehicle (helicopter) having passenger seats 12 
for the purpose of showing how the invention is to be 
used. The invention will display a reading which will be 
indicative of the discomfort experienced by passengers 
in seats 12 during a ride. 

A transducer box 13 containing accelerometers 14-18 
is placed on the vehicle such that accelerometer 14 
produces an electrical signal ejr proportional to the 
vehicle vertical acceleration, accelerometer 15 pro- 
duces an electrical signal ck proportional to lateral ac- 
celeration, accelerometer 16 produces an electrical sig- 
nal e^ proportional to longitudinal acceleration, acceler- 
ometer 17 produces an electrical signal ea proportional 
to roll acceleration, and accelerometer 18 produces an 
electrical signal e^ proportional to pitch acceleration. 
The noise pressure fluctuations interior to vehicle 1 1 
are sensed by a microphone 19 which transforms the 
acoustical fluctuation into an electrical signal e^. 
Electrical signals e;^, 'ey, e„, E^, E^j, and e^y/ are applied ^5 
to a processing unit 20 which is disclosed in detail in 

FIGS. 6(a) and 6(W. 

The electrical signal, cm corresponding to the acous- 
tical pressure fluctuations is applied to the input of an 
octave band analyzer and filter 21 in FIG. 6(a) which ^ 
outputs in digital form the A- weighted noise levels 
within 6 octave bands having center frequencies of 63, 
125, 250, 500, 1000, and 2000 Hz. The six A-weighted 
noise levels, (L^)63, (L^)i25» (L^)250. (L^)soo, (L^)iooo, 
(L/<) 2000 » are then applied as input to a noise discomfort 45 
computer 22 which computes the subjective discom- 
fort, Dm due to the interior noise environment. The 
noise discomfort computer is programmed in accor- 
dance with equations (17) and (18) above to produce 
Dm Note that the noise discomfort computer requires 50 
an additional input, namely the subjective discomfort 
level, D vihi produced by the vibration environment. 

The electrical signals (ew', ey, Cw, ee, and e^,) corre- 
sponding to each of the five vibration acceleration lev- 
els are applied to the input of a fast Fourier transform 55 
signal analyzer 23 which performs spectrum analyses of 
each of the signals and digitally outputs the power spec- 
tral density characteristics (Sty, Sy, Su, S^, S 4 >) of each 
signal. The power spectral density characteristics of 
each vibration are then applied to a spectrum identifica- 60 
tion computer 24 which identifies the random and dis- 
crete frequency components within each spectrum. 
Discrete frequency (or sinusoidal) components are iden- 
tified in terms of the frequency (F) and root-mean - 
square (rms) level {gws, gys» g0j) of each component. 65 
Random components are identified in terms of band- 


by performing a numerical integration of the power 
spectrum between the lower and upper frequency limits 
of the defined bandwidth. The above procedure (steps 
(1) through (4) are repeated for additional peaks lying 
outside of the previously defined bandwidth. Band- 
widths less than 2.0 Hz are treated as sinusoidal vibra- 
tion and band widths equal to or greater than 2.0 Hz are 
treated as random vibration. Additional peaks that lie 
more than 20 dB below the maximum peak are ignored. 
The procedure described above serves to identify the 
random and discrete frequency characteristics of each 
power spectrum. 

The random and discrete frequency characteristics of 
the vibrations are then applied to the appropriate com- 
fort computers 25-31 which transform the physical 
vibration characteristics into subjective discomfort 
units T>wRy Dyj?, D^j?, D„/?, D,^/? for random vibration 
and Dfyv, Dy^, D^j for sinusoidal or discrete frequency 
vibrations. 

Comfort computer 25 is programmed in accordance 
with equation (2) above to produce Dpy/?; comfort com- 
puter 26 is programmed in accordance with equation (3) 
to produce Dy;?; comfort computer 27 is programmed 
in accordance with equation (4) to produce T^br; com- 
fort computer 28 is programmed in accordance with 
equations (la) and (lb) with Table 1(a) to produce D^y^; 
comfort computer 29 is programmed in accordance 
with equations (la) and (lb) with Table 1(b) to produce 
Dys; comfort computer 30 is programmed in accor- 
dance with equations (la) and (lb) with Table 1(c) to 
produce D$s; comfort computer 31 is programmed in 
accordance with equation (6) to produce D^/?; and com- 
fort computer 32 is programmed in accordance with 
equation (5) to produce 

The vertical discomfort values DfVR and Dfvs (out- 
puts of comfort computers 25 and 28) are applied to a 
multiple frequency computer 33 which outputs the 
subjective discomfort, D^y, due to multiple frequencies 
of vibration in the vertical axis. Similarly, Dy/? and Dy^ 
(outputs of comfort computers 26 and 29) are applied to 
multiple frequency computer 34 which outputs the 
subjective discomfort, Dv, due to lateral vibration; and 
Dflft and Das (outputs of comfort computers 27 and 34) 
are applied to multiple frequency computer 35 which 
outputs the subjective discomfort, D«, due to roll 
vibration. The three multiple frequency computers are 
programmed in accordance with equation (19) above. 
The values Dh^, Dy, Da, Du/? and D<^/? are applied to the 
inputs of a combined axis computer 36 which outputs 
the total discomfort, D^, due to the combined axis 
vibrations. The combined axis computer 36 is pro- 
grammed in accordance with equations (7) through 
(12) above. In programming computer 36 


width (Bfy/?, By/^, Bg/?, B„/?, B<f>/?), center frequency 
(CFfy/?, CFy/?, CF$r, CF uRj CF^,/?), and root-mean- 
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The combined axis discomfort, Dc, is then applied to the 
duration correction computer 37 which outputs the 
total duration corrected discomfort due to vibration, 
T>vib’ The duration correction computer is programmed 
in accordance with equations (13) through (15) above. 

The output of the duration correction computer 
is applied to the noise discomfort computer 22 to be 
used in the generation of the noise discomfort contribu- 
tion, Da^. The noise component of discomfort, Dat, is 
then algebrically summed with the vibration compo- 
nent of discomfort, Dy/^, by a summer 38 to generate the 
total subjective discomfort, T>totah attributable to the 
measured noise and vibration environment. The T> total 
signal is then applied to a visual selectable display 39 
which can, in addition to displaying T>totah also displays 
selected values of discomfort produced at earlier stages 
in the meter, i.e., the output display is selectable. Op- 
tions for display include: T> total, ^Vib, ^ noise, ^w, ^v, 
De , t>ur, 

The octave band analyzer and filter 21 which in- 
cludes an analog-to-digital converter, the fast Fourier 
transform analyser 23 which includes an analog-to-digi- 
tal converter, the summer 38 and the selectable display 
39 are all well known and commercially available, and 
hence are not disclosed in detail in this application. 

The principal advantages of this invention over the 
prior art are: Its output is in terms of subjective discom- 
fort units which directly relate to passenger subjective 
acceptance of the measured environment; its output is 
derived from direct software implementation of the 
detailed psychophysical laws governing human subjec- 
tive comfort response to combined noise and vibration; 
its complete recording, analysis, and transformation of 
measured data into subjective units is accomplished on 
board the vehicle and in real time; it is very sensitive to 
slight changes and/or differences in the physical char- 
acteristics (e.g., frequency content, level) of the envi- 
ronment under measurement; it is very useful in design 
tradeoff studies and for comparative evaluations of 
differing ride environments; its effects of multiple fre- 
quencies and multiple axes of vibration are automati- 
cally accounted for within the meter; and it automati- 
cally accounts for the combined effects of noise and 
vibration. 

What is claimed is: 

1. A meter for measuring the ride quality of a passen- 
ger vehicle comprising: 

means on said vehicle for producing an electrical 
signal proportional to the noise level on said vehi- 
cle; 

means on said vehicle for producing several electrical 
signals proportional to vibrations on said vehicle; 
and 

means for combining said electrical signal propor- 
tional to noise and said several electrical signals 
proportional to vibrations to obtain a total discom- 
fort value indicative of the ride quality of said 
passenger vehicle. 

2. A meter for measuring the ride quality of a passen- 
ger vehicle according to claim 1 wherein said means for 
producing several electrical signals proportional to 


vibrations on said vehicle includes several accelerome- 
ters with each producing an electrical signal propor- 
tional to acceleration relative to a single axis. 

3. A meter for measuring the ride quality of a passen- 
5 ger vehicle according to claim 2 wherein said means for 

combining said electrical signal proportional to noise 
and said several electrical signals proportional to vibra- 
tions comprises means receiving said electrical signals 
proportional to vibrations for producing a discomfort 
10 value for each of said accelerations relative to a single 
axis; means receiving said discomfort values relative to 
a single axis for producing a combined axis discomfort 
value; means receiving said electrical signal propor- 
tional to noise level and said combined axis discomfort 
value for producing a noise discomfort value; and 
means for combining said combined axis discomfort 
value and said noise discomfort value to obtain said 
total discomfort value. 

4. A meter for measuring the ride quality of a passen- 
ger vehicle according to claim 3 wherein said means for 
producing a combined axis discomfort value includes 
means for producing a combined axes discomfort value 
that is corrected for time duration. 

5. A meter for measuring ride quality of a passenger 
vehicle according to claim 4 wherein said means for 
producing a noise discomfort value comprises an octave 
band analyzer and filter for dividing said electrical sig- 
nal proportional to noise level into several octave bands 
and a noise discomfort computer receiving said several 
octave bands and said combined axis discomfort value 
that is corrected for time duration for generating said 
noise discomfort value. 

6. A meter for measuring ride quality of a passenger 
35 vehicle according to claim 5 wherein said means for 

producing a discomfort value for each of said accelera- 
tions relative to a single axis comprises a fast Fourier 
transform signal analyzer receiving said several signals 
proportional to vibrations for producing the power 
40 spectral density characteristics of each of said several 
signals; a spectrum identification computer receiving 
said power spectral density characteristics for produc- 
ing values that identify the random and sinusoidal fre- 
quency components within each spectrum; and means 
45 receiving the last mentioned values for producing said 
discomfort value for each of said accelerations. 

7. A meter for measuring ride quality of a passenger 
vehicle according to claim 6 wherein said means for 
producing said discomfort value for each of said accel- 

50 erations includes riieans receiving said values that iden- 
tify the random and sinusoidal frequency components 
within each spectrum for producing said discomfort 
values for each of said accelerations relative to a single 
axis. 

55 8. A meter for measuring ride quality of a passenger 

vehicle according to claim 7 including means for select- 
ably displaying the values: total discomfort value, dis- 
comfort values for the single axis, noise discomfort 
value; combined axis discomfort value; and combined 
60 axis discomfort value corrected for time duration. 

* « * * * 
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